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We have developed a model for neutrino-induced coherent pion production off nuclei in the energy regime 
of interest for present and forthcoming neutrino oscillation experiments. It is based on a microscopic model for 
pion production off the nucleon that, besides the dominant A pole contribution, takes into account the effect of 
background terms required by chiral symmetry. Moreover, the model uses a reduced nucleon-to-A resonance 
axial coupling, which leads to coherent pion production cross sections around a factor two smaller than most 
of the previous theoretical estimates. In the coherent production, the main nuclear effects, namely medium 
corrections on the A propagator and the final pion distortion, are included. We have improved on previous 
similar models by taking into account the nucleon motion and employing a more sophisticated optical potential. 
As found in previous calculations the modification of the A self-energy inside the nuclear medium strongly 
reduces the cross section, while the final pion distortion mainly shifts the peak position to lower pion energies. 
The angular distribution profiles are not much affected by nuclear effects. Nucleon motion increases the cross 
section by ~ 15% at neutrino energies of 650 MeV, while Coulomb effects on charged pions are estimated to 
be small. Finally, we discuss at length the deficiencies of the Rein-Sehgal pion coherent production model for 
neutrino energies below 2 GeV, and in particular for the MiniBooNE and T2K experiments. We also predict flux 
averaged cross sections for these two latter experiments and K2K. 

PACS numbers: 25.30.Pt,13.15.+g,12.15.-y,12.39.Fe 



I. INTRODUCTION 

Neutrino-induced one-pion production off nucleons and nuclei in the intermediate energy region is a source of relevant data 
on hadronic structure. Pions are mainly produced through resonance excitation and these reactions can be used to extract 
information on nucleon-to-resonance axial transition form factors. Besides, a proper understanding of these processes is very 
important in the analysis of neutrino oscillation experiments. For instance, 7r° production by neutral currents (NC) is the 
most important i^,,-induced background to experiments that measure — > ly^ oscillations in the neutrino energy range around 
1 GeV [1]. This is because NC 7r° events can mimic ly^ signal events when, for example, one of the two photons associated with 
the 7r° 77 decay is not detected. This can happen when a photon exits the detector before showering or does not have enough 
energy to initiate a shower. Similarly, 7r+ production by charged currents (CC) is an important source of background in i/^ 
disappearance searches |2l]. 

In reactions on nuclei, pions can be produced incoherently or coherently. In the latter case the nucleus remains in its ground 
state. Coherent reactions are controlled by the nucleus form factor and are more forward peaked than incoherent ones. CC 
coherent pion production has been studied at higher energies in a number of experiments Li 4, 5, 6, 7, 8]. The results could be 
satisfactorily explained by the Rein-Sehgal model which is based on the partially conserved axial current (PCAC) hypothe- 
sis 1 10]. The K2K Collaboration has recently conducted a search for CC coherent pion production induced by muon neutrinos 
with a mean beam energy of 1.3 GeV ifTHl . Contrary to expectations, they found no evidence for CC coherent pion production, 
setting an upper limit of 0.60% for the coherent to total CC pion production ratio. The data show a deficit of forward muons 
in the kinematical region where a sizable coherent production is expected. An attempt to explain this deficit has been done by 
Rein and Sehgal in Ref. 1 12] by including in their model the usually neglected finite muon mass effect ifTsllTil] . In this way they 
find a 25% suppression caused by the destructive interference between the axial vector and pseudoscalar (pion-pole) amplitudes, 
reducing in this way the discrepancy between theory and experiment, though it still persists. This correction affects only CC 
processes and its relevance is reduced as the neutrino energy increases lfT2ll . The negative K2K results are consistent with a very 
recent search performed by the SciBooNE Collaboration | ISfl. 

NC coherent pion production was observed by the Aachen-Padova group lfT6ll on a ^^Al target with both the muon neutrino 
and antineutrino CERN PS beam with average energy of 2 GeV. Positive evidence was also seen by the PS-Gargamelle neutrino 
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and antineutrino Freon experiments lUTtl. V ery recently, the MiniBooNE Collaboration announced the first observation of NC 
coherent n*^ production below 2 GeV flSll. When integrated over the MiniBooNE flux, they find a ratio of coherent plus 
diffractive production over all exclusive NC tt'^ production given by 19.5 ± l.l{stat.) ± 2.5{sys.)% for a mineral oil target 
(CH2). By using Monte Carlo they estimate the coherent rate for a pure ^^C target to be 20.3 ± 2.8{stat.)%. 

On the theoretical side the Rein-Sehgal mod el |9i1 mentioned above assumes that coherent pion production is dominated by 
the divergence of the axial current lfl9[l20l I21I l22ll and can thus be related to the pion-nucleus coherent scattering via PCAC. 
Extrapolation to non-forward angles is done by including a propagator term (1 + Q^/m with niA ~ 1 GeV. The effects 
on the model of considering a finite muon mass was recently analysed in Refs. lfT2ill4ll and, as stated above, they give rise to a 
25% reduction of the CC coherent pion production by muon neutrinos at low neutrino energies. However, one should note that 
the Rein-Sehgal model does not account for nuclear pion absorption, since it does not consider two body mechanisms which are 
those responsible for the absorption of the outgoing pions, and it does not correctly treat quasielastic collisions either Besides, 
the corrections to the outgoing pion angular dependence predicted by the model become quite important for the low neutrino 
energies relevant in MiniBooNE and T2K experiments, as we will show in Subsect. FlVBI 

The PCAC approach was also used in the models of Refs. lEsi I24ll25ll26ll27|] . In Ref. llZTll the authors take into account the 
muon mass effect and include a small non-PCAC transverse current contribution. In all cases the distortion of the final pion 
was included. There are other approaches that do not rely on PCAC. In Ref. [28J coherent pions are produced by virtual A-h 
excitations in the nucleus. The model includes the modifications of the nucleon and A propagators in the medium, evaluated in a 
relativistic mean field approximation, but no final pion distortion was taken into account. KeUcar et al. fl^ improve on the above 
calculation by doing a more sophisticated evaluation of the A self-energy in the medium and treating the final pion distortion 
in a realistic way by solving the Klein-Gordon (KG) equation for a pion-nucleus optical potential. The model of Ref. ll30ll 
uses similar medium corrections and improves on the description of the elementary reaction. On the other hand the final pion 
distortion is treated in the eikonal approximation which is known to fail at low pion energies. In Refs. lf3ll[32ll the authors follow 
the Kelkar et al. calculation in the treatment of the final pion distortion while using, as in Ref. Isoll . a more complete and fully 
relativistic elementary amplitude. 

In a recent publication we have developed a model for CC and NC neutrino- and antineutrino-induced pion production off 
the nucleon in the intermediate energy region [33], which represents the natural extension of that developed in Ref. fsT] for the 
electron analogue eA^ e'N'n reaction. Most previous studies [35. 36. 37. 38,, , 39. 40. 41. 42...43,] of these processes considered 
only the dominant A pole mechanism in which the neutrino excites a A(1232) resonance that subsequently decays into iVvr. In 
our model we have also included background terms required by chiral symmetry (See Fig.[T]below). Some background terms 
were considered before in the works of Refs. [44, 45, 46] although none of these models was consistent with chiral counting. 
In Ref. 1I33I1 . we found that background terms produced significant effects in all channels. As a result we had to re adjust the 
strength of the dominant A pole contribution. The least known ingredients of the model are the axial nucleon-to-A transition 
form factors, of which gives the largest contribution. This strongly suggested the readjustment of that form factor to the 
experimental data, which we did by fitting the flux-averaged v^p jj^pir^ ANL -differential cross section for pion-nucleon 
invariant masses W < lA GeV 1471 148II . Our full model, thus obtained, lead to an overall better description of the data for 
different CC and NC, neutrino- and antineutrino-induced, one-pion production reactions off the nucleon. This reduction of the 
C^(0) value is consistent with recent results in lattice QCD [49], quark model [50] and phenomenological studies [51[. 

Here, we shall apply our model to evaluate CC and NC coherent pion production in nuclei, including the chiral background 
terms in the elementary amplitude. We follow a scheme similar to that advocated in Refs. ll3lll32ll '. but improving the results of 
these latter references by properly taking into account the motion of the nucleons and correcting for some numerical inaccuracies 
that affected the calculations of these two references [52[. Our model should work better close to threshold and hence we 
will concentrate in the neutrino energy range of MiniBooNE and the future T2K experiment where the neutrino peak energy is 
expected to be around 0.6-^0.7 GeV [53]. This work is organized as follows: in sectionllllwe discuss our model for the evaluation 
of CC coherent pion production, including the most relevant aspects of medium corrections for the A and the evaluation of the 
final pion distortion. In section |III] we find the corresponding expressions for the NC coherent pion production case. Finally, in 
sectionHVlwe present and discuss our results. 

II. CC NEUTRINO AND ANTINEUTRINO INDUCED REACTIONS 

We will focus on the coherent CC pion production reaction induced by neutrinos, 

Mk) + Az\g,{pA) ^ r (fc') + Az\gs{PA) + 7r+(fc.) (1) 



' We use a realistic description of tlie A self-energy in the medium and treat quantum mechanically the final pion distortion. 
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The process consists of a weak pion (7r+) production followed by the strong distortion of the pion in its way out of the nucleus. 
In the coherent production the nucleus is left in its ground state by contrast with the incoherent production where the nucleus is 
either broken or left in some excited state. 

The unpolarized differential cross section, with respect to the outgoing lepton and pion kinematical variables, is given in the 
Laboratory (LAB) frame by 



dVL{k')dE'dVt{K) 



Ifcl 4^2%- "^00.+ 



(2) 



with k and k' the LAB lepton momenta, E' = (fc'^ + m^^i/z (jjg energy and the mass of the outgoing lepton, 

G = 1.1664 X 10^^^ MeV^2 the Fermi constant, and E-j^ — [k^ + )i/2 the LAB momentum and energy of the outgoing 
pion, and L and W the leptonic and hadronic tensors, respectively. The leptonic tensor is given by (in our convention, we take 

£0123 = +1 and the metric g^"" = (+,-,-,-)): 



L^fi'J — k^kcr -\- k^k^ — .g^fffc • k + ie^aapk "k^ 
and it is not orthogonal to the transferred four momentum q — k — k' even for massless neutrinos, i.e, L^i'Jq'^ 



(3) 
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FIG. 1 : Model of Ref. jssll for the N —* N'n reaction. The circle in the diagrams stands for the weak transition vertex. 



The hadronic tensor includes all the nuclear effects and it can be approximated by 



64,^3 M2 



-4^+ (g, k,) = j d^fe^'^-'')-'' {ppir) [jp't^ (f; g, k,)] + p,,{f) [j^^^ (f; q, k,)] } 



(4) 
(5) 

with AI the nucleon mass, Pp(„) the nuclear proton (neutron) density, normalized to the number of protons (neutrons). Since 
we have neglected the recoil energy of the final nucleus we have ~ k^{= ^tt)- Finally j7^^+(r; q, k.^), stands for the 
nucleon helicity averaged W^N Ntt^ amplitude evaluated inside the nuclear medium as explained below. Our model for 
the coherent nuclear process is built up from the coherent scattering of the boson with each of the nucleons of the nucleus 
producing an outgoing 7r+. The nucleon state (wave function) remains unchanged in the dispersion and thus after summing 
over all nucleons, we obtain the nuclear densities which appeared in the hadronic tensor of Eq. (|4]l. In the elementary 

W^N Ntt^ process, energy conservation is accomplished by imposing = E^^, while the transferred momentum q — k^r 
has to be accommodated by the nucleon wave functions. Thus, the coherent pion production process is sensitive to the Fourier 
transform of the nuclear density for momentum q — k^r (see Eq. (|5]l). This nuclear form factor gets its maximum value when q 
and are parallel, but for this particular kinematics the vector contribution of the J^j^^ amplitudes, which is purely transverse 
fc^ X q, vanishes. This is the reason why for electron and photon induced reactions, the coherent pion production cross section 
turned out to be a quite small fraction of the total inclusive nuclear absorption one llsillssll . For neutrino induced reactions, the 
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axial contribution of the amplitudes is not suppressed for kinematics where q and fc^r are almost parallel. Thus, the reduction 
induced by the nuclear form factor is much less important, and one might expect a larger relative contribution of the coherent 
pion production channel, as it is the case for some purely hadron reactions (f.i., coherent pion production in the ('^He,t) in 
nuclei |56]). This dominance of the axial contributions has been extensively exploited, through the PCAC hypothesis, to relate 
the neutrino coherent pion production cross section with the pion-nucleus elastic differential one [9,. 12. ,27], 

For the elementary process we have used the model recently derived in Ref. ifssll . In addition to the A pole (AP) mechanism 
(weak excitation of the A(1232) resonance and its subsequent decay into Nn), the model also includes background terms 
required by chiral symmetry. It consists of seven diagrams (see Fig.[T]i: Direct and crossed A(1232)— (first row) and nucleon 
(second row) pole terms (AP, CAP, NP, CNP) contact (CP) and pion pole (PP) contribution (third row) and finally the 
pion-in-flight (PP) term. It provides a fairly good description of all available data for pion production off the nucleon at 
intermediate energies, driven by CC and NC and induced by both neutrino and antineutrino [33.1 . To compute J^^+ (r; g, fc^), 
we will need to evaluate 

i Ur{p')K,N.^Mp)^ * = CAP, NP, CNP, CT, PP, PF (6) 

r 

where the m's are Dkac spinors for the nucleons, normalized such that uu ~ 2M, and the four-vector matrices r^^^+ can be 
read from the explicit expressions of the pion production amplitudes {NTT^\j^^^{0)\N) — J^i ^{p'y^'i-N-n+'^^'P ) ™ ^1- ^-'^^ 
Ref. ll33ll . Finally, p and p' = p + q — k-„ stre the initial and final three momenta of the nucleon. Those momenta are not well 
defined and we approximate the four-momentum of the nucleon N which collides with the by 



p'^ = I JA'P + qY, ] , (7) 



4V " V ' 2 

Hence we assume that the initial nucleon momentum is (/c^ — q)/2 and the final one is — (/ctt ~ <f)/2, with both nucleons being 
on-shell. The momentum transfer is equally shared between the initial and final nucleon momenta. This prescription was firstly 
used in Refs. I54l l55[l for coherent t:'^ photo- and electroproduction, respectively. The approximation is based on the fact that, 
for Gaussian nuclear wave functions, it leads to an exact treatment of the terms linear in momentum of the elementary amplitude. 



In Ref. 1 54] it was shown that this prescription provided similar results as the explicit sum for the nucleon momenta performed 
in Ref. |57]. More recently it has also been employed in Refs. ifssl l59ll for coherent vr" photo- and electroproduction and in a 
recent work on neutrino coherent pion production ll3Tl[32ll . Setting p = —p' = (fc^ — 9 )/2, with energies — given by 
Eq. (|7]i, eliminates some non-local contributions, and it greatly simplifies the sum over all nucleons, which can be cast in terms 
of the neutron and proton densities (see Eq. (|5]l). Furthermore, the sum over helicities in Eq. (|6]l can be also easily performed for 
p = —p' since Ur{p' — —p) — j^Ur{p), so that 



i J2 Mp = -p K,n.+ Mp) = ^Tr [iP + MhX;N.+ ) ' ^ = ^P, CAP , NP . 



(8) 



Thus finally, we will use 



JL+ i^; q,k^)^X. ^v,N.^ i^; q,K), AP, CAP, NP, CNP, CT, PP, PF (9) 

i 

J!;j,^Ar;q,K) = iTr((j6 + Af)7"r^^,+ ) ^, (10) 

expressions that we shall evaluate numerically. Note that the PP term does not contribute to the process since the trace above is 
zero in this case. Within this approximation, the averaged W~^N Ntt^ amplitude inside the nuclear medium, J'^^+ {r; q, fc^r), 
does not depend on r. Below, we will include further medium corrections to the dominant AP mechanism which will induce an 
explicit r dependence. 

Given the importance of the A— pole contribution and since the A properties are strongly modified inside the nuclear 
medium [60, 61, 62, 63], we consider some additional nuclear corrections to this contribution to include the effect of the self- 
energy of the A in the medium S]a(/5(?^))- Here we follow the same approach as in Ref. 113 lH . which is based on the findings of 
Refs. llllllllll . Thus in the A-propagator, we make the substitutions Ma Ma + RbEa and rA/2 r^^"'72 - IhiSa 
and take T,A{p{r)) and r^^""/2 as explained in Sect. II-B of Ref. jlH]. 

So far the formalism has used the bound wave functions of the nucleus, which appear via the proton and neutron densities, and 
has considered only a plane wave for the pion. Pion distortion effects are important, specially for |A:7r | < 0-5 GeV ll29l[30l[3Tl[32ll . 
and are considered here by replacing in Eq. (|5]l 

e-''^---^^(pUir-:k^) (11) 
4e''*^"-^^^iV^:+(f;4) (12) 
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The pion wave function </?* + (r; k-^) corresponds to an incoming solution of the Klein Gordon equation, 

Ifil+ir; k^) = El^l+{r; k^) , (13) 

with Vopt (r) the optical potential which describes the 7r+ -nucleus interaction. This potential has been developed microscopically 
and it is explained in detail in Refs. [63, 64]. It contains the ordinary lowest order optical potential pieces constructed from the 
s- and p-wave ttN amplitudes. In addition second order terms in both s- and p-waves, responsible for pion absorption, are 
also considered. Standard corrections, as second-order Pauli re-scattering term, angular transform term (ATT), Lorentz-Lorenz 
effect and long and short range nuclear correlations, are also taken into account. This theoretical potential reproduces fairly well 
the data of pionic atoms (binding energies and strong abso rptio n widths) [631 and low energy 7r-nucleus scattering ||64l[ . At low 
pion energies, it is an improvement over the one used in ll3lll32ll . that was based on A dominance of the vrA^ interaction. Another 
possible improvement would be the inclusion of the Coulomb interaction between the outgoing pion and the nucleus. This can 
be taken into account by means of the replacement 

E^^E^-Vcir^ (14) 

in the right hand side of Eg. dT3]) . where Vc {r) is the Coulomb potential created by the nucleus, including finite size and vacuum 
polarization effects, see ['6J, o?[. We will discuss the effect of this correction in Sect.lIVI 

The replacement in Eq. ( fTSl i. that takes into account the fact that the pion three momentum is only well defined asymptotically 
when the pion-nucleus potential vanishes, induces some non-localities in the amplitudes. To treat these non-localities we have 
adopted the following scheme: 

• In Eq. (|4l), we approximate {k^^l, which arises from the phase space integrations, by the modulus of the asymptotic three 
momentum (i?^ — m^)2 

• In the AP, CAP, NP, CNP terms, we note that there exist either a NNn or a NAtt vertex (see Eq. (51) of Ref. iH), 
which induces a factor fc" in the amplitudes. Indeed, for those terms we could re-write 

J!;N^Ar;q,K)^iK)aJ!:^^Ar;q,k^), i ^ AP,CAP ,NP ,CNP . (15) 

We do not consider any non-locality in the tensor J7^^^^+, and we use the prescription of Eqs. ( fTTT l and ( fT2] l to account 

for /ctt in the contraction between k" and j'^^^^ in Eq. ( fTSl ). This approach to treat the non-localities is equivalent to that 
assumed in refs. ll3ll[32ll . 

• We do not consider any non-locality for the CT and PP contributions. 

Antineutrinos induce the coherent production of negatively charged pions. To study these processes, we use ||33|[ 

L^;} = 4;^ (16) 

'7p';-[„.-](^;9,fc^) = J,'^^+[p^+]{r;q,k^) (17) 

and implement the appropriate changes in the pion-nucleus Vopt and Vc potentials to properly account for the distortion of the 
outgoing 7r~ ll64l[ . 

Differences between neutrino and antineutrino induced cross sections are proportional to the interferences among the axial 
and vector current contributions. Since the latter ones are suppressed by the nuclear form factor, as we discussed after Eq. (|5]), 
we expect roughly similar neutrino and antineutrino cross sections. This will also be the case for the NC driven processes studied 
in the next section. 



-V + ml + 2E^Vopt{r) 



III. NC NEUTRINO AND ANTINEUTRINO INDUCED REACTIONS 

To extend the above formulae to the case of NC 7r° coherent production, 

Mk) + AzlgsiPa) ^ Mk') + Azlgsip'a) + 7r°(fc.) (18) 

we have. 



dn{k')dE'dn{K) ~ \k 1 16^2 Nc.o 



(19) 



(g, fc.) = / e'«--~ {pp(f) [j-^^o (f; q, k^] + pjr) [j^^^, (f; q, k^] } (f; K) (21) 



6 



with k' and E' = \k' \ the LAB outgoing neutrino momentum and energy. The leptonic tensor is given in Eq. (O and it is now 

orthogonal to = (fc — k')'^ for massless neutrinos, i.e, L\i^}q^ = L^^Jq"^ = 0. 

Both lepton and hadron tensors are independent of the neutrino lepton family, and therefore the cross section for the reaction 
of Eq. (dill is the same for electron, muon or tau incident neutrinos. Furthermore, the hadron tensor is the same for neutrino and 
antineutrino induced reactions, and thus to study antineutrino reactions we just have to change the sign of the antisymmetric part 
of the leptonic tensor (see Eq. (fTSIl). 

To evaluate the hadronic tensor, we use the model for the NC pion production off the nucleon derived in Ref. ifssll and thus 
we have 

Ji^^a (f ; g, fc^ ) = ^ Ji:^^o (r ; q,k^), i = AP, CAP, NP, CNP, CT, PP, PF (22) 

i 



(f; q,k^) = -Tr[{P + M^' r^^^„ j - . (23) 

with the same prescription for the nucleon momentum as in the CC case. Within this model the PP, PF and CT diagrams do 
not contribute to the NC ir^ production off the nucleon. The AP and CAP terms provide equal Z^p — > pvr^ and Z'^n rni^ 
amplitudes, with F^p.^y^,, and T^^p.^^o obtained from T^p.^^^ and Tp^p p^^ multiplying these latter matrices by the 

overall factors 2\/2/(3 cos 9c) and 2\/2/ cos 9c respectively, and multiplying the vector form factors by (1 — 2 sin^ 9w), being 
9c the Cabibbo angle and 9w the Weinberg angle. Direct and crossed nucleon pole terms lead to 

pA" _ _i nA^P „ p + i + M f^^^ ^^^^ T^NP_( 1 Ptt" 



^ CNP-Ntt° 



■,D^^^I±{v,%iq)-A-,.^iq) 



— Klb, D 



CNP 



with ~ 93 MeV the pion weak decay constant, = 1-26 the axial nucleon coupling, and 



1 ptt" 



— 1 nTT 







Z:N 



2 X 



Ft{q')l^+l 



4" — 



2M 



-cr qu 



N 



Giiq')ri5 + (Giiq') + G\{q'))^^ + Gf,(r j ) ^75 



(25) 

(26) 
(27) 



N 



The pseudoscalar part of the axial current, which is proportional to q^^, does not contribute to the differential cross section for 
massless neutrinos. Besides the Z^NN form factors are given by [.65i1 



ifizF^Y''" = ±MvPj'-2sin2 0M/Mp,„Pf"--^,P2^ 



(C^) ' = ±Ga — 



(28) 

(29) 
(30) 



For isoscalar nuclei, the direct and crossed nucleon pole terms do not contribute because the existing cancellation between neu- 
tron and proton contributions, and we have total dominance of the A mechanisms. If we neglect the vector current contributions^, 
finite lepton mass effects and approximating cos 9c ~ 1, we find that the CC coherent pion production cross section is twice 
the NC one, as deduced from the relevant isospin factors and the factor of four of difference between Eqs. (|2|i and (fT9] l, for CC 
and NC driven processes, respectively^. For non symmetric nuclei, as long as the A dominance holds, we will reach the same 
conclusion. Nevertheless, we remind here that for low and intermediate muon neutrino energies (< 1.5 — 2 GeV), one should 
expect sizable corrections (25% at 1.3 GeV, and greater at smaller energies) to the approximate relation crcc ~ 2ctnc due to the 
finite muon mass 1IT2I1 . 

To evaluate pion distortion effects we compute the vr" — wave function by using the appropriate pion-nucleus optical poten- 
tial ll64ll and setting Vc = 0. Non-localities in the amplitudes are treated as in the CC case. 



^ They will be suppressed by the nuclear form factor (see discussion after Eq. (5))- 

^ For the case of isoscalar nuclei, the approximate relation acc ~ "^(^NC is far more general and it can be directly deduced from PCAC and isospin invaiiance. 
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Nucleus 


Tp [fm] 


r„[fm] 


a 




1.692 


1.692 


1.082 


16 Q 


1.833 


1.833 


1.544 


20«Pb 


6.624 


6.890 


0.549 fm 



TABLE I: Charge (rp, a) and neutron matter (Vn, a) density parameters for different nuclei as given in Ref fe^l. For carbon and oxygen we use a modified 
hai'monic oscillator density, p(r) = PQ{l + a{r/rt,i)'^) exp(— (r/rjv)^), while for lead, we use a two-parameter Fermi distribution, p(r) = po/(l + exp((r — 
rjv)/a)). 



IV. RESULTS 

We shall always use the full model of Ref. ll33ll where the dominant nucleon-to-A axial form factor was fitted to data 
resulting in Crf (0) = 0.867 and Ma\ = 0.985 GeV. Note that the Goldberger-Treiman relation, traditionally assumed in the 
literature, implies a larger value of C^(0) 1.2. We will come back to this point below, at the end of the discussion of the 
results shown in Fig.|4] 

Firstly, we compile in Table Uthe input charge densities, taken from Ref. fbd\, used in this work. For each nucleus we take the 
neutron matter density approximately equal (but normalized to the number of neutrons) to the charge density, though we consider 
small changes, inspired by Hartree-Fock calculations with the density-matrix expansion ll67ll and corroborated by pionic atom 
data [68]. However, charge (neutron) matter densities do not correspond to proton (neutron) point-like densities because of the 
finite size of the nucleon. This is taken into account by following the procedure outlined in section 2 of Ref. [68.1 (see Eqs. 
(12)-(14) of this reference). 



A. General Results 

First, in Fig.|2]we show the pion momentum distribution (LAB) for CC and NC coherent pion production induced by and 

on a ^*^0 (CC case) and ^^C (NC case) targets. In the upper panels we show the CC case for a z^^, beam energy of 600 MeV, 
which is in the expected peak energy region of the future T2K experiment. In the lower two panels we show NC results for ?ivv 
beam energy of 850 MeV. In all panels, the short-dashed line corresponds to our results in plane wave impulse approximation 
(PWIA), in which we use a plane wave for the outgoing pion and neglect medium effects on the dominant A contribution. 
Medium effects introduced through the A self-energy play a very important role largely reducing the PWIA results. This is 
shown by the long-dashed line in the upper-left panel. When the pion distortion is also taken into account (via the substitutions 
of Eqs. (fTTT i and ( fT2] i) the cross section is further reduced, and the peak is shifted towards lower energies reflecting the strong 
absorption and the higher probability of a quasielastic collision of the outgoing pion by the nucleus in the A kinematical region. 
The total cross section reduction is around 60% for a beam energy of 600 MeV. Our full model results thus obtained are shown 
by the solid line. Medium and pion distortion effects in coherent pion production were already evaluated in Refs. ||3TI [32ll . 
However in these works the motion of the nucleon was neglected. This is to say, though the authors of these references also 
use the prescription p = —p' = (fc^ — <f)/2 to compute the elementary W^{Z^)N nN' amplitude, however the nucleon 
momenta in the Dirac's spinors appearing in Eq. (jSJ are neglected. The effect of putting the nucleons at rest can be clearly seen 
in the dotted line of the upper-left panel and results in a 15% decrease of the total cross section. Results obtained for the 
antineutrino CC induced reaction are shown in the upper-right panel. The cross section is some 30% smaller than in the neutrino 
case due to the sign change in the axial part of the lepton tensor which results in a different vector-axial interference. Similar 
effects are seen for NC reactions shown in the lower panels of Fig. |2] In this case the antineutrino cross section is reduced by 
just 10% with respect to the neutrino one. 

In Fig. |3]we show the effect of Coulomb distortion on the outgoing charged pion. One can see the expected shift in the peak 
towards higher (lower) energies of the positive (negative) pion distribution when the Coulomb distortion is taken into account. 
The net effect in the total cross section is nevertheless small, amounting to a 5% change for beam energies in the 500 MeV 
region. For higher energies the effect is expected to be less important. 

We now study, for CC reactions, the effect of including the background terms on top of the dominant direct (AP) and crossed 
(CAP) A pole contributions. As mentioned before, the PF term does not contribute to the coherent cross section. Besides, the 
direct (NP) and crossed (CNP) nucleon pole term contributions partially cancel each other, while the chiral background terms 
CT, PP vanish for isospin symmetric nuclei due to an exact cancellation between proton and neutron contributions. This latter 
cancellation is partial for asymmetric nuclei. In the left panel of Fig. |4]we show the results for ^^C. As seen in the figure the 
effect of the background terms, both in the PWIA and in the full calculation (including medium effects and pion distortion), are 
very small, thus corroborating the findings of Ref. [32]. In the right panel of the figure we show full calculation results for ^osp^ 
, which is the most asymmetric nucleus with possible experimental interest. In this latter case the inclusion of the background 
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FIG. 2: Pion momentum differential cross section in the LAB frame for different coherent pion production reactions. Short-dashed line (in blue) has been 
calculated using planes waves for the outgoing pion and without including any in-medium correction for the A. Results with A nuclear medium effects ai'e 
shown in the upper-left panel by the long-dashed line (in green). Our full model calculation, including medium effects on the A and the distortion of the outgoing 
pion wave function, is shown by the soUd line (in red). Finally, the effect of putting the nucleons at rest is shown in the upper-left panel by the dotted line (in 
magenta). 



terms reduces the cross section in an appreciable way. We find similar conclusions for NC driven processes. 

This predominant role played by the A mechanism, in conjunction with the findings of Ref. [33] and the fact that the coherent 
pion production reaction in nuclei is mostly driven by the axial part of the interaction, allows us to conclude that most of the 
previous theoretical studies \3V\ of the pion coherent channel might be overestimating the cross section roughly 

by a factor of two. This can be easily understood as follows. Background terms turned out to be very important at the nucleon 
level and because of them, the flux-averaged i/^p fj^^pn^ ANL ^^-differential cross section i47i |48I1 is described with an 
axial form factor C^(0) of around 0.9 133], significantly smaller than the traditionally used value of about 1.2, deduced from 
the Goldberger-Treiman relation. This reduction of the contribution of the A pole mechanism in the weak pion production 
off the nucleon is compensated by the non-resonant terms. However, when one studies the neutrino coherent pion production 
in isoscalar nuclei we find a negligible contribution of the non-resonant terms, and thus the cross section is determined by the 
axial part of the A mechanism of which gives the largest contribution. Thus, we predict cross sections around a factor of 
(1.2/0.9)^ ~ 2 smaller than those approaches which assume the Goldberger-Treiman relation"^. This fact was firstly pointed 
out by the authors of Ref. [32], who used for the very first time the background terms derived in Ref. [33] for the elementary 
reaction. However, we improve here the results of this reference by properly taking into account the motion of the nucleons, as 
discussed above, and correcting for some numerical inaccuracies [52] (of the order of 20% in the total cross section, and larger 
at the peak of the da/dk^r differential distribution) that affected the calculations of this reference and those of a previous work 



A word of caution is required here. There exists some degree of inconsistency among the ANL (A^JaS] and BNL f6^ measurements of the integrated u^p — ► 
/i~p7r+ cross section, being the latter larger than the former one. The model of Ref. |33], including non-resonant background terms, with C^(0) = 1.2 
(consistent with the off diagonal Goldberger-Treiman relation) would lead to a better description of the BNL data than if the lower vale of around 0.9 is used. 
Thus, if one favours BNL data one could still use a high value of (0) = 1.2, which would lead to larger coherent pion production cross sections. 
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FIG. 3: Coulomb distortion effects on the pion momentum differential LAB cross section for CC coherent pion reactions induced by Ufj,, i/^. 




FIG. 4: Pion momentum differential LAB cross section, calculated with and without background terms, for CC coherent pion production induced by a 500 
MeV Ufj, beam. The curves labeled as "Only A" stand for results obtained from the AP and CAP mechanisms. Left panel: results for a ^^C target and PWIA 
and full calculations; Right panel: results for a ^"*Pb target. 

Next we pay attention to the differential cross section, and in Fig.|5]we show this distribution for CC driven processes in 
the energy region of the future T2K experiment. There, we also show da/dq^^^, where q^^^ is calculated, under the assumption of 
a Quasi-Elastic (QE) neutrino-nucleon interaction, from the measured outgoing lepton energy and scattering angle in the LAB 
frame, 

ql,^~2Mq^^^ - -2M {El'^^ - E') 
^ ME' - 

M - E' + \k'\cos9' 

The q'^^^ distribution as compared with the q'^ one is clearly shifted to lower absolute values and peaks roughly at 0. This fact 
might be used to reduce the QE background by requiring that coherent events should have a reconstructed q^^^ value smaller 
than some appropriate cut, as was done in the K2K analysis carried out in Ref. [11]. 

We turn now to angular distributions of the final pion and muon in induced CC reactions. In Fig. |6]we show the pion 
angular LAB distribution with respect to the incoming neutrino direction. As expected, and due to the nucleus form factor, the 
reaction is very forward peaked. Inclusion of medium effects on the A propagator and the final pion distortion largely reduce 
the cross section. The angular distribution profile keeps its forward peaked behaviour, although less pronounced once the pion 
distortion is included. This can be seen on the right panel where the angular distributions are all equally normalized to one. Such 
a behaviour can be understood by taking into account that the pion wave function (^* (r; k^r) has not well defined momentum, 
in contrast with the plane wave e"'*^" ''. Putting the nucleons at rest has some effect on the cross section but hardly affects the 
angular distribution profile. The situation is very similar for the muon angular distribution, shown in Fig.|2]for CC coherent pion 
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production by a 0.65 GeV energy beam on a ^^C target. In this case the angular distribution profile is completely unaffected 
by the nuclear effects on the A and the pion distortion. 



B. The Rein-Sehgal Model and the NC MiniBooNE E^{1 - cos 6^ ) distribution 



We have also examined the NC differential cross section with respect to the variable Et^{1 — cos 6*^), proposed by the Mini- 
BooNE Collaboration in its recent analysis of coherent tt*^ production of Ref. [18]. Pion variables are referred to the LAB frame. 
Results are shown in the left top panel of Fig. [8] for a neutrino beam energy of 800 MeV (close to the energy peak of the 
MiniBooNE experiment) on carbon. Our differential cross section is appreciably narrower than that displayed in Fig. 3b of 
Ref. |18]^. Indeed, our distribution at — cos^tt) = 0.1 GeV has already fallen off by a factor 18, while the MiniBooNE 
distribution has fallen off by less than a factor 4 at 0.1 GeV, and even at 0.2 GeV the reduction factor is still smaller than 15. 
We find hard to understand this discrepancy. Since the excitation of A(1232)— resonance mechanism is dominant, we expect 
the outgoing pion to have a total energy of around 0.25 — 0.35 GeV (see for instance left bottom panel of Fig.|2l). On the other 



^ Note however, that the i?.,r(l — cos 0^) distributions shown in Ref. 1 18] are not proper differential cross sections. This is because, they have not been corrected 
for acceptance or cut efficiencies and are plotted for reconstructed kinematic quantities. So they include the effects of the selection criterion (the efficiency of 
which can vary as a function of E^{1 — cos 9^)), as well as reconstruction effects in the MiniBooNE detector. Currently, the MiniBooNE Collaboration is 
working to have all of the effects of the detector, event reconstruction, and selection removed. The new results will follow in an upcoming paper, where actual 
differential cross sections will be available L70i1 . 
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FIG. 7: Muon angular differential cross section for CC coherent pion production induced by a 650 MeV energy i/^ beam on a ^^C target. Left panel: absolute 
differential cross section; Right panel: differential cross sections normalized to one. The 6p angle is refeiTed to the incoming neutrino direction in the LAB 
frame. 



hand, we find that the cross section is almost negligible for cos > 0.7, as can be inferred from the pion angular LAB differ- 
ential cross section displayed in the right top panel of Fig.|8] Thus, we easily understand why we find that the E^ril — cos6'7r) 
distribution becomes quite small above 0.1 GeV. To have non-negligible signals in the E.^{1 — cos 6'^) — 0.1 — 0.2 GeV region 
would require, approximately, values of cos 6*^ in the interval 2/3 — 1/3, which would translate into values of 9-^ in the range 
50° — 70°. Since g is strongly aligned to the incoming neutrino direction, such high pion angles look hardly compatible with the 
forward character of the coherent reaction, which is just imposed by the nucleus form factor (Fourier transform of the nuclear 
density for momentum g — fc^r) and the q ■ k.^ dependence of the amplitudes (dominated by the axial contribution). 

The strong disagreement between our prediction and the MiniBooNE histogram in the region below £'^(1 — cos6'7r) < 0.1 
GeV is even most worrying, because there, the cross sections are much larger We are aware that within the MiniBooNE 
flux, there exist neutrino energy components higher and smaller than the 800 MeV considered in the top panels of Fig. [8] In 
the left bottom panel of Fig. |8] we show the E-,^{\ — cos^jr) differential cross section for three more neutrino energies (300, 
550 and 1300 MeV), in addition to that of 800 MeV considered in the top panels. In all cases, we find very small signals for 
i?7r(l — cos (?^) above 0.05 GeV. Also in this panel, we show our £'7^(1 — cos 6-^) differential cross section convoluted with the 

MiniBooNE flux (solid line), and we certainly find a distribution definitely narrower than that published in Ref. [18[]. 

Since the MiniBooNE analysis relies on the Rein-Sehgal model for coherent vr" production H, the strong shape difference 
should be understood in terms of the differences between our model and that of Ref. 19D. 



I. The t~dependence of the Rein-Sehgal model 



Rein and Sehgal made use of the Adler's PCAC formula OlOtl and approximated (both for neutrino or antineutrino induced 
processes) the coherent tt" production differential cross section by 



(33) 



with X = —q^/2Mq^, y = (f jEi,, and t = {q — k-^)'^ = —{q— A:^)^. Besides, the nuclear form factor is calculated as 

{pp{r) + p„(r )}, and finally Fabs is a <— independent attenuation factor^. The above expression was 
deduced in the so called parallel configuration, for which the fc^ and k'^ four momenta are proportional (therefore q^ — 0) and 

cos 9qk^ (angle formed by fc^ and q) and |fc^|/|g | are approximated to one everywhere except in the nuclear form factor. It was 
continued to non-zero q^ values by including a propagator term of the form (1 — q^/m^j)^^, with niA ~ 1 GeV. The model 



* In the original work of Ref. it is stated that -Fabs takes into account effects of pion absorption in the nucleus. As defined in Ref. K3], ^abs only removes 
from the flux pions that undergo inelastic collisions but, as explained below, no true absorption is actually included. 
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should work well close to this parallel kinematics, and constitutes a good approximation at the high neutrino energies, above 2 
GeV, explored in the original work of Ref. |9]. However, the approximations in which is based this model become less justified 
as the neutrino energy decreases. For the energies relevant in the MiniBooNE and T2K experiments, non parallel configurations 
turn out to be more important, and the Rein-Sehgal model predictions are less reliable. 

Actually, we see from Eq. (|33] | that the Rein-Sehgal differential cross section depends on cos 0^ or t only through the nuclear 
form factor, and any further cos and/or t behaviour induced by the dependence of the amplitudes on k-^ is totally neglected. 
However, it is reasonable to expect that these additional angular dependences might play a role when the pion emission is not 
completely forward. To illustrate this point, we have re-derived Eq. (l33T l from our model. We have considered the dominant 
axial contribution of the AP mechanism and made use, in this case, of the Goldberger-Treiman relation to express C^(0) in 



terms of the ttA^A coupling, /*, and the pion decay constant /^r C^{0) = 



Besides, we have considered neither 



the in medium A— selfenergy, nor the pion distortion effects. On the other hand, we have used the TVAtt Lagrangian^ 

C^NA = ^^^f\d^$)^ + h.c. (34) 

to compute da{'K^N — > 7r°Af)/d|t|. To further simplify, we have worked in the non-relativistic limit for the baryons. Apart 
from -Fabs, we reproduce Eq. ( l33T l. but with an extra (cos Oqk^ Y factor in the right hand side of the equation, which for the case 
of parallel kinematics is equivalent to (cos OttY . At high neutrino energies, the pion is emitted strongly forward and thus, it is 
consistent to approximate this factor by one, as it was done in the original work of Ref. |9] to obtain Eq. ( [33] ). For MiniBooNE 
neutrino energies, relatively large 6*^ angles are allowed, in special for low pion momenta ( | /ctt | < 0.2 GeV) [71], and this factor 
(cos6't)^ could make the E.^{\ — cos 6*^) distribution significantly narrower than that predicted by the Rein-Sehgal model. 
For a fixed value of Et^{1 — cosO.,^), the effect becomes more important as the pion energy decreases. For instance, if we fix 
£'^(1 — cos^tt) = 0.05 GeV, we find that (cos^tt)^ ~ 1/2 or 3/4, for an averaged pion energy of 0.17 or 0.375 GeV, respectively. 
As it was shown in the talk of J. Link at Nuint07 JtH, in the Rein-Sehgal model, low energy pions < 0.5 GeV) produce 
wider ( 1 — cos Oj^ ) shapes than those of higher energies ( | fc^ | > 0.5 GeV), and thus we expect this type of corrections to be 
important. In addition, there also exist corrections that vanish in the ^ Q (or equivalently \q\/q^ ^ 1 ) and/or [fcTrl/lff | ^ 1 
(also implied by the t = Q approximation in the amplitudes) limits. 

Another way to see the limitations of the Rein-Sehgal model is the following. As mentioned, this model assumes no further 
dependence on t than that encoded in the nuclear form factor. Since t = implies q = fc^, we have replaced k" in Eq. ( fTSl l 
by q". It is to say, we replace k^ by q in the pion emission vertex**. To better compare with the Rein-Sehgal predictions, we 
have again just considered the dominant axial contribution of the AP mechanism, without considering pion distortion and 
A in medium effects, MiniBooNE flux convoluted results are displayed in the right bottom panel of Fig. |8] We see that 
the new E-^{1 — cos6'7r) distribution is significantly wider than that obtained without implementing this replacement, and that it 
reasonably describes the MiniBooNE published distribution (solid histogram in the plot). The agreement is much better, when 
we compare with some preliminary MiniBooNE results (dashed histogram) obtained with a different treatment of the outgoing 
pion Final State Interaction (FSl), as we will explain in the next sub-subsection^. Without giving a special meaning to this 
agreement, this simple calculation serves the purpose of illustrating the uncertainties associated to the t — approximation at 
low energies, for which the nuclear form factor still allows some deviations from the completely forward scattering. 

We conclude that, the Rein-Sehgal pion coherent production model for MiniBooNE and T2K experiments is not as reliable as 
for the case of neutrino energies above 2 GeV. We expect sizable corrections to the predictions of this model, both for differential 
distributions and for integrated cross sections. Our model provides an i?7r(l — cos 9.^) distribution much more peaked around 
zero, and thus it might improve the description of the first bin value in Fig. 3b of Ref. [IS]. Moreover, the drastic change in 
the i?7r(l — cos^Tr) distribution shape might produce some mismatch between the absolute normalization of the background, 
coherent and incoherent yields in the MiniBooNE analysis. One the other hand, and besides of the issue of the used value'° 
for C^(0), the Rein-Sehgal model, adopted by the MiniBooNE Collaboration in Ref. fil8l1 . overestimates by a large factor the 
coherent integrated cross section for MiniBooNE energies. This is due first to the t = approximation in the amplitudes assumed 
in this model that produces a too wide E-^{1 — cos 6,^) distribution, which leads to cross sections larger by about a factor of two 
than those obtained when the t— dependence is properly taken into account (see for instance the different areas below the solid 



^ Here, (/> is the pion field, 4'^ is a Rarita Scliwinger = 3/2"*" field, is tlie isospin transition operator (vector under isospin rotations and its Wigner-Eckart 

irreducible matrix element is taken to be one) from isospin 1/2 to 3/2. 
* If we were to repeat with this replacement the derivation of Eq. (33) as explained above, we will recover it exactly (apart from -F^bs), and the correction 

factor (cos fqfe^ )^ would not appear. 
' We are indebted with G. Zeller for providing us these preliminary results. 

" Within the Rein-Sehgal model, this constant is implicitly fixed to approximately 1.2, since the Goldberger-Treiman relation is used to express the coherent 
production cross section in terms of the elastic tt^^N n^^N one. 
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and dashed-dotted curves in the right bottom panel of Fig. [8]l. Secondly, because all sort of in-nuclear-medium effects, like pion 
absorption or modification of the elementary ttN ttN cross section" inside of the nucleus..., which are not accounted for in 
the Rein-Sehgal model, and that turn out to be relevant at MiniBooNE energies. 

2. The {1 - cos 0^) > 0.1 GeV region and NUANCE FSI 

We turn now to the region i?^ ( 1 — cos 6'^ ) > 0. 1 Ge V. Also in this region, the bulk of the discrepancies among our predictions 
and the MiniBooNE results can be understood in terms of the inaccuracies of the Rein-Sehgal model at low energies, as can 
be appreciated in the right bottom panel of Fig. [8] But here, the outgoing tt'^ FSI effects, incorporated by the analysis of 
Ref. [18], induce some additional discrepancies. The MiniBooNE analysis relies on the Rein-Sehgal models for incoherent 
(resonant) [72] and coherent [9] tt*^ production, which are implemented in the NUANCE event generator [73]. In the case of 
coherent production, the Rein-Sehgal model includes an absorption factor Fabs to effectively account for the pion wave function 
distortion from a plane wave (see Eq. ([33])). This Glauber-type factor decreases exponentially with the pion-nucleon inelastic 
cross section (see Eq. (24) of Ref. [9]) and thus it neither accounts for pion absorption, which is a two nucleon mechanism, 
nor it does for quasielastic distortion. Quasielastic steps, induced by elastic pion-nucleon collisions and not forbidden by Pauli 
blocking, excite and/or break the nucleus, and are not removed by this factor Fabs- This Glauber factor only removes events 
where the 7r° suffers an inelastic collision with a nucleon, and as a result changes its charge, or other mesons are produced in the 
final state. 

The procedure followed in Ref. ifTsIl to describe coherent pion production is somewhat different. They set Fabs to 1 in the 
Rein-Sehgal model of Ref. f^, and implement absorption as part of the FSI. In our understanding, coherent pion production 
cross section cannot be calculated from a Monte Carlo cascade algorithm. This is because, by definition, the coherent production 
is a one step process'^, and the quantum mechanical transition matrix element gives the amplitude probability for producing a 
pion outside of the nucleus, which is left unchanged. The coherent contribution should be incoherently added to that due to the 
inelastic channels to find out the total pion production cross section. 

Nevertheless, one could still reasonably estimate the total coherent cross section from the NUANCE FSI cascade, if it would 
be used to eliminate from the flux of outgoing neutral pions, not only those which get absorbed or those that suffer inelastic 
processes, like multiple pion production, meson production or pion charge exchange, but also those that undergo quasielastic 
steps, induced by elastic pion-nucleon collisions, in their way out of the nucleus. However, to our knowledge, these latter 
events are accounted for in the MiniBooNE analysis, despite they are not coherent since the final nucleus, as a result of the 
secondary coUisions, is not left in the ground state. We believe this is acknowledged by the authors of Ref. ifisll when they say 
"...that rescattered events with a tt'^ in the final state may be misclassified in NUANCE, as would be the case when a coherently 
produced 7r° rescatters elastically through a resonance" . As a result of these collisions, the 7r°'s might change their direction 
and give rise to events in the NUANCE cascade at significantly larger values of 0^. Indeed in the right bottom panel of Fig. [8] we 
observe significantly less events above F^(l — cos 6'^) > 0.1 GeV when the NUANCE FSI is turned off (dashed histogram)'^. 
The effects below 0. 1 GeV are much smaller, and in total, the change in the shape leads to a reduction of around 20% in the 
integrated cross section. 

In our calculation of the coherent cross section, we certainly remove those secondary events by means of the optical potential 
employed to compute the pion wave function. The imaginary part of the pion-nucleus potential is responsible for the removal 
of flux of the outgoing pions on their way out of the nucleus. This imaginary part is due to pion absorption, but also to pion 
quasielastic steps. Hence, the use of the full optical potential will eliminate the pions which are absorbed and also those which 
scatter quasielastically. We might try to theoretically estimate this effect by switching off the quasielastic contribution to the 
pion-nucleus optical potential induced by elastic pion-nucleon collisions, and using an optical potential with an imaginary part 
due to absorption and inelastic channels alone. In this way, we will remove the absorbed pions and those that undergo inelastic 
colhsions, but not those which scatter quasielastically, which will still go out of the nucleus and are accounted for by the 



' Within our model, we include these modifications by means of the consideration of the A— selfenergy, by considering the Pauli blocking in the computation 
of A— decay width in the medium, and by using the pion wave function, ^* , instead of a plane wave. 

^ The nomenclature here might be confusing. There exist multiple step contributions to the coherent reaction. For instance, a A is formed in a NC scattering, 
it decays with the nucleon falling back into the hole created by A formation, the decay -k^ creates a subsequent A, which in turn decays emitting a tt" that 
escapes the nucleus and the associated nucleon also drops into the ground state configuration. The point we want to make here is that such contributions 
cannot be taken into account in a Monte Carlo cascade algorithm. This is because it would require the coherent sum of the multiple step amplitudes, while a 
Monte Carlo alogorithm uses probabilities (cross sections). We do include these multiple step contributions within our formalism thanks to the use of a pion 
wave function solution of the Klein-Gordon equation, with an optical tt" -nucleus potential, instead of using a plane wave. 

^ Note that when NUANCE FSI is turned off, besides of getting rid of the unwanted quasielastic steps, pion absorption is not taken into account. However, this 
latter effect though, produces a diminution of events, it does not significantly change the shape of the Et^{1 — cos ) distribution. On the other hand, since 
the histogram has been re-scaled down, the overall normalization is not an issue any more. Nevertheless, we should point out that, there could be same minor 
differences in the acceptance or cut efficiencies with respect those used in the published histogram fT3l . 
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MiniBooNE Monte Carlo generator. This was considered in the past in the study of the pionic decay of A— hypernuclei I*74[ p75ll . 
leading to moderate enhancements of the decay widths of the order of 10 — 15% in ^^C |i76jl . We find here similar effects, and 
for the MiniBooNE flux averaged cross section, we find an enhancement of around 20% (see Table HHi. in good agreement with 
the effects observed by turning off the NUANCE FSl. 



C. K2K, MiniBooNE and T2K Flux Averaged Cross Sections 

In Table HIl we show our predictions for the K2K ifTlll and MiniBooNE ifTsIl flux averaged cross sections as well as for the 
future T2K experiment. In Fig.|9] we show some results for K2K and MiniBooNE experiments. In all cases, we normalize the 
neutrino or antineutrino flux (j) to one. In principle, we would like to compute the corresponding convolution with the neutrino 
or antineutrino flux 

relish 

dE(j)'{E)a{E), i = K2K, MiniBooNE, T2K (35) 

with iJj^Q^, iJJjjgjj the lower and upper flux limits, and a{E) the corresponding CC/NC muon/electron neutrino/antineutrino 
induced nuclear coherent cross section, as a function of the neutrino/antineutrino energy^'*. In practice, the predictions of our 
model become less reliable when the energy increases, since the model neglects all resonances above the A(1232). Sophisticated 
recent calculations, as those of Refs. |30, 31, 32], suffer from exactly the same limitation. That is the reason why we have set up 
a maximum neutrino energy {E^^^) in the convolution, and approximated 



J^hdEcj,^{E)a{E) 
J^HdE^iE) 



(36) 



where we fix the upper limit in the integration (neglecting the long tail of the neutrino fluxes) to i?max = 1-45 and 1.34 GeV for 
CC and NC muon neutrino/antineutrino driven processes, respectively. The phase space for the fifth differential 



dn{k')dE'dn{k^) 

cross section, up to irrelevant constants, is determined by |fc'||fc7r| [see Eqs Q, (EJl and ( fT9l ), (|20tl. For CC processes, for 
instance, and muon neutrino energies of around 1.45 GeV, the phase space peaks at pion energies of around 730 MeV, which 
leads, neglecting the nucleon momentum, to vrA^ invariant masses below 1.5 GeV. Up to these energies, one can reasonably 
assume A(1232) dominance. 

In the case of the K2K experiment a threshold of 450 MeV for muon momentum was imposed as an additional selection 
criterion flT]. We have implemented this cut also here, and in that case we have been able to go up to £'j„^^^^^=1.8 GeV. In 
these circumstances, we still cover about 90% of the total flux in most of the cases. For the T2K antineutrino flux, we cover just 
about 65% of the total spectrum, and therefore our results for the convoluted cross sections are less reliable. 

For neutrino energies above 1 GeV, and though the A contribution plays a central role in pion production fl^, one should 
bear in mind that other resonances could certainly be also important. This would affect the results presented in Table and 
Fig lD Taking into account that the T2K and MiniBooNE fluxes peak at neutrino energies of around 0.6-0.7 GeV, where the A 
resonance contribution is much more dominant, is reasonable to expect corrections (higher cross sections) of around 20-30% to 
our results for these two experiments. Certainly, the corrections could be larger for the K2K case, since it that case the neutrino 
energy spectrum peaks at higher energies, around 1.2 GeV. 

We see that our prediction, subject to some uncertainties, lies well below the K2K upper bound, mainly thanks to the use of 
a low value for C^(0), while we predict a NC MiniBooNE cross section notably smaller than that given in the PhD thesis of 
J.L. Raaf [77]. However, this latter value should be taken with extreme caution. It was obtained from a preliminary analysis, 
that since then has been notably improved. Moreover, the MiniBooNE Collaboration has not given an official value for the 
total coherent cross section yet, and only the ratio coherent/(coherent+incoherent) has been presented 1 18]. Nevertheless, as we 
have discussed at length, we believe the MiniBooNE analysis might overestimate this ratio, not only because some of the tt^'s 
which undergo FSI collisions are accounted for as coherent events instead of being removed, but more importantly because the 
Rein-Sehgal model predicts an incorrect (wider) E^^il — cosOtt) shape for coherent 7r"'s. The first of the effects produces an 
enhancement of the coherent cross section of the order of 20% (see the MiniBooNE NC* entry in the table), while it is much more 
difficult to quantify the second of the effects. This is because, in addition to the variation of the integrated area, it might produce 



Note that the cross section trivially vanishes for neutrino/antineutrino energies below the pion production threshold, which obviously is different for CC and 
NC driven processes because of the final lepton mass. 
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a possible mismatch between the absolute normalization of the background, coherent and incoherent yields in the MiniBooNE 
analysis. To finish this discussion, we would like to point out that the K2K cross section and the value quoted in Ref. |77] seems 
somehow incompatible with the approximate relation crcc ~ SctnC; expected from A— dominance and neglecting finite muon 
mass effects (see discussion at the end of Sect.HIIll. 

Our predictions are about 20-30% smaller than those obtained in Ref. ll32ll from model II, which uses our value C^(0)« 0.9. 
This discrepancy is due to different facts: the existence of some numerical inaccuracies affecting the results of Refs. lIsTi [32I1 . 
the inclusion of the non zero momenta in the nucleon spinors, different ranges in the flux convolution, etc... 

For the future T2K experiment, we get cross sections of the order 2—2.6 x lO^^'^cm^ in carbon and about 2.3— 3.0 x 10^*°cm^ 
in oxygen, considerably smaller than those predicted by the Rein-Sehgal model. 

The SciBooNE Collaboration has recently set 90% confidence level upper limits on the cross section ratio of CC coherent 
pion production to the total CC cross section at 0.67 x 10~^ at mean neutrino energy 1.1 GeV and 1.36 x 10~^ at mean neutrino 
energy 2.2 GeV II15I1 . If we use a value of 1.05 x 10"'^^ cm^/nucleon for the total CC cross section, as quoted in [15], the 
SciBooNE upper limit for the ratio will transform in a upper bound of about 8 x lO^'"' cm^ for the coherent cross section 
at 1.1 GeV. At 1.1 GeV, our prediction for the CC coherent pion production cross section in carbon is 5.7 x 10~*° cm^ (see 
points in the left panel of Fig.|9]l, which is totally compatible with the SciBooNE bound. However, according to the Rein-Sehgal 
model JH [T2l implemented in the SciBooNE simulation, the cross section ratio of CC coherent pion production to the total CC 
cross section is calculated to be 2.04 x 10^^ in Ref. ifisll . The SciBooNE limits correspond to 33% and 67% of the Rein-Sehgal 
model prediction at 1.1 GeV and 2.2 GeV, respectively f\3]. From our discussion in Subsect. [IV BK see last paragraph in the 
sub-subsection [TV B It , we easily understand why the Rein-Sehgal model overestimates the coherent cross sections by a large 
factor at neutrino energies around 1 GeV (approximately a factor of two because of the value'^ of C^(0), and approximately 
another factor of two because of the t = approximation in the amplitudes, in addition of all sort of in-nuclear-medium effects, 
like pion absorption or modification of the elementary ttN ttN cross section in the medium,... which are not accounted for in 
the Rein-Sehgal model). We also understand why at the higher neutrino energy of 2.2 GeV the Rein-Sehgal model works better, 
since the larger the energy, the better the t — approximation in the amplitudes and smaller the nuclear effects become. Note 
also, that at 2.2 GeV, we expect heavier resonances than the A(1232) to play an important role, and thus the issue of the value 
of C^(0) is less relevant. 

To conclude this section, in Fig.[TO]we show muon neutrino/antineutrino CC and NC coherent pion production off carbon 
and oxygen targets. We see that both for CC and NC driven processes, the ratio of neutrino over antineutrino cross sections 
approaches one as the neutrino energy increases. This is due to the nuclear form factor which reduces the vector contribution to 
the amplitudes (and therefore the interference between the vector and the axial parts) as the neutrino energy increases. Besides, 
and as a consequence of the nuclear form factor and other nuclear in medium effects, we also see that cross sections do not scale 
as A^, being A the nuclear mass number, as expected from a coherent reaction. This is in good agreement with the findings of 
Refs. [9, 31-[. Indeed at 1 GeV, oxygen and carbon cross sections turn out to be in a proportion of around 6 to 5, instead of 1.8 
to 1, as it would be deduced from an A^— type scaling law. Finally, we observe sizable corrections to the approximate relation 
CTcc ~ 2tTNC for these two isoscalar nuclei in the whole range of neutrino/antineutrino energies examined in this work. As 
pointed out in Refs. lfT2l [Till , this is greatly due to the finite muon mass, and thus the deviations are dramatic at low neutrino 
energies. In any case, these corrections can not account for the apparent incompatibility among the CC K2K cross section and 
the NC value quoted in Ref. iIttIi . mentioned above. 

V. CONCLUSIONS 

We have developed a model for neutrino/antineutrino CC and NC coherent pion production off nuclei which is based on a 
microscopic model for neutrino/antineutrino-induced one-pion production off the nucleon derived in Ref. [33]. The model of 
Ref. [33[ includes the dominant A production mechanism, but it also takes into account background terms required by chiral 
symmetry. In its application to coherent production we have further taken into account the main nuclear effects expected to be 
important in reactions off nuclei. While the model presented here is similar to the one in Refs. [^3Tl l32|] . we have improved on 
that calculation by taking consistently into account the nucleon motion, and by using a more sophisticated pion optical potential. 
The consideration of the nucleon motion increases the cross section by a non negligible amount, while Coulomb effects on 
the emission of charged pions lead to small changes in the cross section. Moreover, we have corrected for some numerical 
inaccuracies ]'5?] (of the order of 20% in the total cross section, and larger at the peak of the da/dky^ differential distribution) 
that affected the calculations carried out in Refs. ll3Tll32ll . 

In agreement with Refs. 1I3II I32I1 . we find a strong reduction of the cross section, mainly due to the modification of the A 
self-energy in the nuclear medium, and a shift to lower energies of the outgoing pion distribution due to the final pion distortion. 



Note, however that, (0) will partially cancel out in the ratios measured by the SciBooNE Collaboration. 
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-cose,) (GeV) cos6». 




FIG. 8: Laboratory i?^ (1 — cos 9,) and cos 9-,^ distributions for the u^^C u^^Cn^ reaction, at MiniBooNE energies. Different models are considered 
in the upper panels, while our full model is always used in the left bottom one, where we also show the i?,r(l — cos 6,) distribution convoluted with the 
MiniBooNE flux (solid line). Details of the convolution are explained in the text (see Eq. )36t and TablellTl. In the right bottom panel, we show results from the 
axial contribution of the AP mechanism, neglecting pion distortion and A in medium effects (see the text for the explanation of the two curves). In both 
bottom panels, we display the MiniBooNE published histogram (solid), conveniently scaled down, taken from the right panel of Fig. 3 in Ref. 1 18]. Finally, in 
the right bottom panel, we also show MiniBooNE results (dashed histogram) obtained by turning off the NUANCE ESI of the outgoing pion (G. Zeller private 
communication) . 



The angular distributions of both pion and muons with respect to the incoming i/^ direction are forward peaked due to the nuclear 
form factor While the muon angular distribution profile is almost unaffected by nuclear corrections, in the pion case, part of 
the strength is shifted to larger angles due to the distortion of the final pion wave function. Non-resonant terms, that turned out 
to be very important at the nucleon level llssll '^. give small contributions to the coherent pion production off isospin symmetric 
nuclei. This leads us to find coherent pion production cross sections around a factor of two smaller than most of those previously 
published. 

We have also performed a detailed discussion of the MiniBooNE results of Ref. Hisll and the analysis performed there to 
identify NC coherent 7r° events. We have shown that the Rein-Sehgal model used in this analysis is not accurate enough in 
this case. This is because the MiniBooNE flux mainly consists of neutrinos below 2 GeV, and for such low neutrino energies, 
the corrections to the outgoing pion angular dependence predicted by the Rein-Sehgal model become quite important. As a 
consequence, the Rein-Sehgal model leads to distributions notably wider and integrated cross sections much larger than those 



Their inclusion made necessary to reduce the nucleon-to-A resonance axial coupling. 
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TABLE II: NC/CC muon neutrino and antineutrino coherent pion production total cross sections for K2K, MiniBooNE and T2K experiments. In the case 
of CC K2K, the experimental threshold for the muon momentum | fc' | > 450 MeV is taken into account. To convert the cross section ratio given in fllll into 
a coherent cross section (K2K), we use the value of 1.07 X 10""^* cm^/nucleon for the total CC cross section, as quoted in 1 11]. For the MiniBooNE NC* 
entry, we present our results when an optical pion-nucleus potential with an imaginary part due to absorption and inelastic channels alone is used to compute 
the distortion of the outgoing pion (see text for more details). The absolute NC -k^ coherent cross section quoted in the PhD thesis of Ref. 1 77] should be taken 
with extreme caution, since in the published paper (Ref. 118]) it is not given. There, it is quoted the ratio of the sum of the NC coherent and diffractive modes 
over all exclusive NC tt" production at MiniBooNE. Some details on the flux convolution are compiled in the last three columns. 



12^ 12^ 

V C -> V C -K 



1 ' 1 ' 1 ' 


1 ' 1 


• O(E^) [10"*°cmT 




— (t)(E^) [GeV"'] 




— o(E^)(])(E^) [lO'^cm^GeV"'] 






T2K 


^ I.I.I. 


^ ^ . 1 — 



0.4 



0.6 



0.8 



1.0 



1.2 



1.4 



(GeV) 



o(E^) [10 cm'] 

(1)(E^) [GeV"'] 

o (E ) (|) (E ) [10"''°cm"GeV"'] 



MiniBoone 




FIG. 9: CC (left) and NC (right) coherent pion production cross sections in carbon. We also show predictions multiplied by the T2K (left) and MiniBooNE 
(right) neuti'ino energy spectra. In the region of neutrino energies around 0.6 GeV, the lower curves stand for the T2K and MiniBooNE Vpi fluxes nonnalized 
to one. 



predicted in this work. Finally, we have predicted muon neutrino/antineutrino CC and NC coherent pion production off carbon 
and oxygen up to neutrino energies of the order of 1.4 GeV, and convoluted those cross sections with the K2K, T2K and 
MiniBooNE fluxes. Our cross sections are considerably smaller than those predicted by the Rein-Sehgal model. 

We expect the present model to provide accurate coherent pion production total and differential cross sections in the first 
resonance region, where the A(1232) plays a relevant role. This energy region is very important for the analysis of present and 
forthcoming neutrino oscillation experiments for which good and reliable theoretical calculations are needed. 
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